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Complete spectra of the staggered Dirac operator jp are determined in four-dimensional SU{2) gauge fields 
with and without dynamical fermions. An attempt is made to relate the performance of multigrid and conjugate 
gradient algorithms for propagators with the distribution of the eigenvalues of jp . 



1. INTRODUCTION 

Big efforts have been undertaken to find effi- 
cient multigrid (MG) methods for the computa- 
tion of propagators in background gauge fields; 
see the lists of references in [1,2], and also [3]. 
Although ultimately one wants to simulate theo- 
ries with dynamical fermions, all previous works 
focussed only on quenched gauge fields. How- 
ever, it is reasonable to expect that MG meth- 
ods have a chance to perform better when one 
considers "real" gauge fields which are gener- 
ated in the presence of dynamical fermions. On 
the other hand one will not expect any big dif- 
ference for the behavior of the conjugate gradi- 
ent (CG) algorithm. The reasons for these two 
statements are as follows. The inclusion of the 
fermionic determinant in the Monte Carlo pro- 
cess will tend to decrease the number of (approx- 
imate) zero modes. This is so because configu- 
rations with less low-lying modes are more prob- 
able. MG methods intend to take care of the 
low-lying modes (which are responsible for criti- 
cal slowing down) on coarser grids, and the task of 
dealing with a reduced number of low-lying modes 
should be easier. Concerning the CG algorithm, 
its (asymptotic) convergence properties are de- 
termined by the condition number. Since condi- 
tion numbers of the (negative squared) massless 
staggered Dirac operator are not influenced dra- 
matically by the presence of dynamical quarks, 
one does not expect a significant consequence for 

*to appear in the proceedings of LATTICE'94 

tWork supported by Deutsche Forschungsgemeinschaft. 



the convergence behavior of CG. In the present 
study we focus on the consequences of dynamical 
fermions for the performance of a variational MG 
algorithm which proved to work in the quenched 
case, but which was unable there to outperform 
CG [4]. 

2. MULTIGRID METHOD 

The variational MG method for the solution of 

(-p 2 + m 2 ) X = f (1) 

is described in detail in Refs. [1,4]. We use a 
blocking procedure which is consistent with the 
symmetries of free staggered fermions. In this 
scheme coarsening is done with a factor of 3 [5]. 
Since the volumes of subsequent layers of the MG 
differ by a factor of 81, the actual implementaion 
of the MG method was only a twogrid algorithm. 
The kernel of the averaging operator is defined 
as the solution of a gauge covariant eigenvalue 
equation. In Ref. [5] our choice was given the 
name "Laplace choice", because we project on a 
block-local approximation of the "fermionic two- 
link Laplacian" A. This A is defined through 
p 2 = A + cr^F^v, where is the lattice defi- 
nition of the field strength by means of plaquette 
terms. 

3. SPECTRUM OF -p 2 

As explained in the introduction, naively one 
expects p to have less approximate zero modes 
in the presence of dynamical fermions than in the 
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quenched case. In order to study this conjec- 
ture we need firstly a Hybrid Monte Carlo pro- 
gram, and secondly a method to determine the 
low-lying spectrum of p. For the generation of 
four-dimensional SU(2) gauge fields coupled to 
dynamical staggered fermions a FORTRAN pro- 
gram with vectorized CRAY code was used, which 
had been written by S. Meyer and B. Pendleton. 
They used this program when they studied the 
chiral transition with many fermion flavors in the 
SU{2) Higgs model [6]. Meyer's and Pendleton's 
program was used with four flavors of staggered 
fermions. The spectrum of p was determined by 
means of a Lanczos procedure. 

3.1. Lanczos procedure 

The Lanczos method has been used in lattice 
field theory for a long time, see e. g. [7]. In the 
present exploratory study the complete spectrum 
of — p 2 was determined. This was done in or- 
der to be sure about the correctness (i. e. to have 
no numerical uncertainties) in the distribution of 
low-lying modes. Of the several computational 
variants of the Lanczos procedure the most sta- 
ble one as described in [8, Algorithm 9.2.1 and 
remark on p. 492] was implemented. The "good" 
eigenvalues were determined by means of Cul- 
lum's and Willoughby's procedure [9]. More de- 
tails can be found in Ref. [10]. 

3.2. Numerical results 

One knows a priori that every eigenvalue of 
— p 2 in SU(2) gauge fields is fourfold degener- 
ate [11]. Thus, if there are no further degenera- 
cies, then on a lattice A of volume |A| there must 
be |A|/2 different eigenvalues. Their sum must 
equal ^Tr (—p 2 ) = 4|A|. These statements are 
valid for periodic and for antiperiodic boundary 
conditions. 

Concrete numerical investigations were done 
with the following parameters. The spectrum of 
— p 2 was investigated on 6 4 and 12 4 lattices. Two 
different kinds of boundary conditions (b.c.) were 
used. One with periodic b.c. in all directions for 
gauge and Fermi fields, and one with periodic b.c. 
for the gauge field in all directions and with an- 
tiperiodic b.c. on the Fermi field in time direc- 
tion and periodic b.c. in spatial directions. The 
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Figure 1. Spectrum of — p 2 on 12 4 lattices with 
antiperiodic boundary conditions for the Fermi 
field in t-direction. 

coupling (3 = 4/gr 2 of the Wilson action for the 
SU(2) gauge fields was varied between 1.8 and 
5.0. Quark masses m in the Hybrid Monte Carlo 
runs were chosen to be m = 0.2 and m = 0.05 (in 
units of a = twice the spacing of the staggered 
lattice). These values were also used in Meyer's 
and Pendleton's work [6]; they quote m = 0.1 and 
m = 0.025 since they measured physical quanti- 
ties in units of a/2. 

It turned out [10] that in nontrivial gauge fields 
there seem to be no additional degeneracies to 
the ones explained above. Only for (3 = 5.0 
it was impossible to identify |A|/2 eigenvalues 
whose sum equals 4|A|. For (3 = 1.8, . . . , 2.8 Cul- 
lum's and Willoughby's method works perfectly. 
On 6 4 (12 4 ) lattices we always found 648 (10368) 
eigenvalues whose sum came out as 5184 + ^6 
(82944 + <5i 2 ), with S 6 < 8- 10" 9 (S 12 < 1.7- 10" 6 ) 
in REAL arithmetic on a CRAY Y-MP. Be- 
cause of the randomness of the nonvanishing off- 
diagonal matrix elements of — p 2 this is very good 
evidence that the spectrum was determined ex- 
actly. 

Spectra on 12 4 lattices are shown in Fig. 1. 
We number the different eigenvalues by k = 
0, 1,2,.. ., with Ao < Ai < . . . The data shown 
for finite (3 are results obtained with a gauge field 
generated by the Hybrid Monte Algorithm in the 
presence of dynamical fermions with a mass of 
m = 0.2. However, on the scale of the whole 
spectrum there is very little difference compared 
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to m = 0.05 or to a quenched gauge field. For 
finite (3 there is also no difference between peri- 
odic and antiperiodic b.c. 1 The (3 = oo values 
are shown in Fig. 1 with their true multiplicities 
modulo the fourfold degeneracy mentioned above, 
so that the (3 = oo values indicate the curve which 
the numerical data should approach for large (3. 

When one looks at the complete spectra, fig- 
ures for 6 4 lattices look practically the same on 
the overall range as Fig. 1, when one rescales the 
abscissa by 16. (Only for the eigenvalues of the 
free — p 2 this is not true.) 

Results for low-lying spectra and for condition 
numbers can be found in Ref. [10]. 

4. INVERSION OF (-p 2 + m 2 ) 

A result of the present study is that the con- 
vergence behavior of CG in nontrivial gauge fields 
is practically only determined by the condition 
number k of {—p 2 + m 2 ), and by the lattice 
size. This should be the case when the eigen- 
values are distrubuted uniformely between the 
lowest and the highest one [12]. For configu- 
rations on a lattice of given size with the same 
ac, CG yields sequences of RMS norms of residu- 
als which practically coincide, even if the spectra 
are different. As mentioned above, on the over- 
all range of the spectra there is little difference 
between quenched simulations and simulations 
with dynamical fermions (of mass m = 0.2, 0.05). 
Therefore slight fluctuations in the distrubution 
of eigenvalues on small scales do not affect the 
convergence of CG. Thus if one wants to study 
the convergence of the CG algorithm one can do 
that with "cheap" quenched gauge fields, one does 
not have to take "expensive" unquenched config- 
urations. 

Results for the two-grid algorithm are as fol- 
lows. An obvious statement is that convergence 
of the MG algorithm is not determined by k. This 
is clear in the limiting case of free fields, because 
in pure gauges critical slowing down is completely 
eliminated by MG, i.e. convergence is completely 
independent of k. In nontrivial gauge fields con- 

1 This should be the case for expectation values of gauge 
invariant quantities when —11 is an element of the gauge 
group. I wish to thank R. Sommer for this remark. 



vergence of MG depends on details of the spec- 
trum. Unfortunately, as in quenched gauge fields, 
inferiority of MG was found compared to CG, a 
factor of about 10 in CPU time. The poor per- 
formance of MG found earlier [4] is no feature 
of quenched computations. More details are in 
Ref. [10]. 
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the CRAY Y-MP of the Univ. of Kaiserslautern. 
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